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ABSTRACT 

Quark clustering could occur in cold quark matter because of the strong coupling 
between quarks at realistic baryon densities of compact stars. Although one may still 
not be able to calculate this conjectured matter from first principles, the inter-cluster 
interaction might be analogized to the interaction between inert molecules. Cold quark 
matter would then crystallize in a solid state if the inter-cluster potential is deep 
enough to trap the clusters in the wells. We apply the Lennard-Jones potential to 
describe the inter-cluster potential, and derive the equations of state, which are stiffcr 
than that derived in conventional models (e.g., MIT bag model). If quark stars are 
composed of Lennard-Jones matter, they could have high maximum masses (> 2M Q ) 
as well as very low masses (< 10 _3 M Q ). These features could be tested by observations. 
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1 INTRODUCTION 

To understand the nature of pulsars we need to know the 
state of cold quark matter, in which the dominant degree 
of freedom is quarks, and their Fermi energy is much larger 
than their thermal energy. However, this is a difficult task 
because of (i) the non-perturbative effect of the strong in- 
teraction between quarks at low energy scale and (ii) the 
many-body problem of vast assemblies of interacting parti- 
cles. 

On one hand, some efforts have been made for under- 
standing the behavior of quantum chromo-dynamics (QCD) 
at high density, among which a color super-conductivity 
(CSC) state is currently focused on in pertu rbative QCD 
as we ll as in QCD-based effective models (e.g.. lAlford et all 
2008). On the other hand, it is phenomenologically conjec- 
tured that astrophy sics! cold quark matter could be in a 
solid state |Xu|[2003), since the strong interaction may ren- 
der quarks grouped in clusters and the ground state of re- 
alistic quark matter might not be th at of Fermi gas (see a 
recent discussion given bv lXull2009l ). If the residua inter- 
action between quark clusters is stronger than their kinetic 
energy, each quark cluster could be trapped in the potential 
well and cold quark matter will be in a solid state. Solid 
quark stars still can not be ruled out in both a strophysics 
and particle physics (|Horvathll2005l : IOwenll2005l ). Addition- 
ally, there is evidence that the interaction between quarks is 
very strong in hot quark-gl uon plasma (i.e ., the strongly cou- 
pled quark-gluon plasma, IShu rvak 2005^), according to the 
recent achievements of relativistic heavy ion collision exper- 
iments. When the temperature goes down, it is reasonable 
to conjecture that the interaction between quarks should be 
stronger than that in the hot quark-gluon plasma. 



Because of the difficulty to obtain a realistic state equa- 
tion of cold quark matter at a few nuclear densities, we 
try to apply some phenomenological models, which would 
have some implications about the properties of QCD at low 
energy scale if the astronomical observations can provide 
us wit h some limitatio ns on such models. In our previous 
paper (|Lai fc Xvj||2009h . a polytropic quark star model has 
been suggested in order to establish a general framework 
in which theoretical quark star models could be tested by 
observations. This model can help us to understand the ob- 
served masses of pulsars and the energy released during some 
extreme bursts; however, this is a phenomenological model 
and does not include the form of interaction between quarks. 
To calculate the interaction between quarks and predict the 
state of matter for quark stars by QCD calculations is a dif- 
ficult task; however, it is still meaningful for us to consider 
some models to explore the properties of quarks at the low 
energy scale. 

We can compare the interaction between quark clus- 
ters with the interaction between inert molecules. A single 
quark cluster inside a quark star is assumed to be color- 
less, just like each molecule in a bulk of inert gas is elec- 
tric neutral. The interaction potential between two inert gas 
molecul es can be well descr ibed by the Lennard-Jones po- 
tential (jLennard- Joneslll924I ) 



u(r) 



4t/o[(^) 12 
r 



(1) 



where Uo is the depth of the potential and ro can be con- 
sidered as the range of interaction. This form of poten- 
tial has the property of short-distance repulsion and long- 
distance attraction. We assume that the interaction between 
the quark clusters in quark stars can also be described by 
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the form of Lennard- Jones potential. Q If the inter-cluster 
potential is deep enough to trap the clusters in the poten- 
tial wells, the quark matter would crystallize and form solid 
quark stars. Under such potential, we can get the equation 
of state for solid quark stars, where the pressure comes from 
both the inter-cluster potential and lattice vibrations. Be- 
cause the chromo-interaction is stronger than the electro- 
magnetic interaction that is responsible for the intermolec- 
ular forces, the values of parameters, Uo and ro used in the 
cold quark matter should be different from that in the inert 
gas, and should be determined in the context of quark stars. 

The model of quark stars composed of Lennard- Jones 
matter is much different from the conventional models (e.g., 
MIT bag model) in which the ground state is of Fermi gas. In 
the former case the quark-clusters are non-relativistic parti- 
cles, whereas in the the latter case quarks are relativistic par- 
ticles. Consequently, the equations of state in this two kinds 
of models are different, and we found that the Lennard- Jones 
model has some more stiffer equations of state, which lead to 
higher maximum masses for quark stars. Certainly, a quark 
star can be very low massive (< 10 _3 Mq) due to self color 
interaction. On the other hand, we find that under some 
reasonable values of parameters, a quark star could also be 
very massive (> 2Mq). 

This paper is arranged as follows. The details of lattice 
thermodynamics are listed in Section 2. The forms of equa- 
tions of state are given in Section 3, including the compari- 
son with the MIT bag model, and we show the corresponding 
mass-radius curves in Section 4. We make conclusions and 
discussions in Section 5. 



2 LATTICE THERMODYNAMICS 

For an inter-cluster potential which is deep enough to trap 
the clusters in the potential wells, the quark matter would 
crystallize to get lower energy and form solid quark stars. 
In this section we use the results in classical solid physics to 
discuss the properties of crystallized cold quark matter. 



2.1 The inter-cluster potential 

Like the inert gas, the interaction potential u between two 
quark-clusters as the function of their distance r is described 
by the Lennard- Jones potential, Eq.(l). Let us consider a 
system containing N clusters with the volume V, then the 
total interaction potential is 



U 



N 



U 



(2) 



and if we ignore the surface tension, we get 



1 The interaction bet ween nuclei can be described by the a — ui 
model jWaleckalll974h , which is also characterized by the short- 
distance repulsion and long-distance attraction. Recently, the 
nucleon -nucleon potenti al has been studied by lattice QCD simu- 
lations lllshii et al.11200'^) , and they also derives a strong repulsive 
core at short distances. The interaction between quark-clusters 
in cold quark matter could also have long-distance attraction 
and short-distance repulsion if scale and vector mesons contribute 
there. We note that this short-distance repulsion is essential to 
reproduce a stiff equation of state in our model. 



{4t/ [(^) 12 -(^) 6 ]}. (3) 



The lattice structure of cold quark matter is unknown, and 
we adopt the simple-cubic structure. The cold quark mat- 
ter may have other kinds of structures, but that will not 
make much differences at least quantitatively. If the near- 
est distance between two quark-clusters is R, then the total 
interaction potential of N quark-clusters is 



U(R) = 2NU [A 12 (^) 12 - A 6 (^)% 



(4) 



where A12 = 6.2, and Aq = 8.4. In the simple cubic struc- 
ture, the number density of clusters n is 

n = R~ 3 , 



(3) 



U(R) = 2NU {A 12 rl 2 n A ~ A 6 rf )n 2 ), (6) 
and the potential energy density is 

e p = 2U (A 12 r 1 2 n 5 - A 6 r 6 n 3 ), (7) 
2.2 Lattice vibrations 

Consider a system of volume V containing TV quark-clusters. 
Each quark-cluster in the crystal lattice undergoes a three- 
dimensional vibration about its lattice site. Performing a 
normal-mode analysis in which the vibrations of the lat- 
tice are decomposed into 3N independent normal modes of 
vibrations, the total lattice vibration is a superposition of 
these 3N decoupled vibrations. 

For cold quark matter, the thermal vibration can be 
neglect compared to the zero-point energy of phonon, so 
the average energy of an individual mode of vibration with 
frequency uij is 

Ei = ±nu s . (8) 

Using Debye approximation, at low temperature, the ther- 
modynamics properties of crystal lattice is mainly deter- 
mined by the long-wavelength sound waves, and the prop- 
agation of the wave can be decomposed into one longitude 
mode with velocity v\\ and two transverse modes with veloc- 
ity v±, and the total velocity v has the relation 



1 , 1 



3 v 



and the total energy of the 3N vibrations is 



(9) 



(10) 



where f(cu)duj is the number of modes in the interval from 
lu to uj + du), and ui m is the maximum frequency related to 
the non-continuous structure of the solid, and is determined 
by 



E = I -fiujf(uj)dw, 



la 2 sl/3 

ujm = v(bir n) 



(11) 



Under the condition of zero-temperature, the integration can 
be done and the total energy of the crystal vibrations is 



— 9V , 2,1 4 
E = — (6tt )3 7wm3 , 
8 

so the energy density of the lattice vibration is 



(12) 
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-(67r 2 )3R w 3. 



(13) 



3 EQUATION OF STATE FOR QUARK STARS 

3.1 Quark stars composed of Lennard-Jones 
matter 

The total energy density for cold quark matter is 

2 



= e P + £l + nm c c 
= 2U (A 12 rl 2 n 5 - A 6 r 6 n 3 ) 

+ — (67T 2 ) 3 hvnz + nm c c 2 , 



(14) 



here m c is the mass of each quark-cluster. The pressure can 
be derived as 



P* = 



dn 



= 4U (2A 12 r 1 2 n 5 ~A e rln 3 ) + ^{6n 2 )3hvn^. (15) 

8 

Apart from quarks, there are electrons in quark matter. 
In MIT bag model, the number of electrons per baryon N B /A 
is found for different strange qu arks mass m s and coupling 
constant a s dFarhi fc Jaffe|[l984h . In their results, when a s — 
0.3, N e /A is less than 10~ 4 ; a larger a B means a smaller N e /A 
at fixed m s , because the interaction between quarks will lead 
to more strange quarks and consequently less electrons. In 
our model, we also consider the strong interaction between 
quarks as well as between quark-clusters, and consequently 
the required number of electrons per baryon to guarantee 
the neutrality should be also be very small. Although at 
the present stage we have not got the exactly value for the 
number density of electrons, we assume that N D /A is less 
than 10 -4 . We find that even at this value, the pressure 
of the degenerate electrons is negligible comparing to the 
pressure of quarks. Therefore, we neglect the contribution 
of electrons to the equation of state. Then the equation of 
state for quark stars is 



P = P« 

P = £q/C , 

3.2 Parameters 



(16) 
(17) 



Up to now, there are a couple of parameters in the equation 
of state, and in this section we will show how to determine 
them. 

1. The long- wavelength sound speed v for lattice vi- 
bration. For the extremely relativistic systems, the sound 
velocity is l/v3, and in general it will be less than l/v3- 
We find that the equation of state does not change much 
when v goes from the velocity of light c to 10~ 5 c, so in our 
calculations we set v = c/3. 

2. The mass of quarks. Quark stars are composed en- 
tirely of deconfined light quarks (up, down and strange 
quarks), the so-called strange stars. Because the deconfined 
phase transition and the chiral-restoration phase transition 
might not occur simultaneously in the QCD phase-diagram, 
we give each quark a constituent mass and assume it is one- 
third of nuclear mass. 

3. The number of quarks in one cluster N q . Quarks 



are fermions and they have three flavor (up, down and 
strange) degrees of freedom and three color degrees of free- 
dom. Pauli's exclusion principle tells us that if in the inner 
space quarks are exchange-asymmetric, they are exchange- 
symmetric in position space and they have a tendency to 
condensate in position space. We therefore conjecture the 
existence of quark-clusters in quark matter and leave the 
number of quarks in one cluster, N a , as a free par ameter. A 
18-quark cluster, called quark-alpha ((Michel 1988), could be 
completely asymmetric in spin, flavor and color space, so in 
our calculation we set A?q=18, and we also set jV q =3. On the 
other hand, it has been conjectured that strongly interacting 
matter at high densities and low temperatures might be in a 
"quarkyonic" state, which also contains three quarks in one 
cluster, a nd is characterized by chiral symmetry and con - 
finement l|McLerran fc Pisarski[|2007l ; iBlaschke et al.ll2008l ). 
However, in our model, the state is characterized by chiral 
symmetry breaking and deconfinement. 

4. The depth of the potential Uq. Given the density of 
quark matter p and the mass of each individual quark, from 
Heisenberg's uncertainty relation we can approximate the 
kinetic energy of one cluster as 



E k 



lMeV(f)§(§; 
Po 18 



(18) 



where po is the nuclear matter density. To get the quarks 
trapped in the potential wells to form lattice structure, Uo 
should be larger than the kinetic energy of quarks. Because 
of the strong interaction between quarks, we adopt (7o=50 
MeV and 100 MeV to do the calculations. 

5. The range of the action ro. Since a quark star could 
be bound not only by gravity but also by strong interaction 
due to the confinement between quarks, the number density 
of quarks on a quark star surface p s is non-zero. For a given 
p s , we can get ro at the surface where the pressure is zero. 
We choose p B as 2 times of nuclear matter densities, and get 
the value of ro accordingly, which is found in the range from 
about 1 to 3 fm. 

When Uo and ro are given, the inter-cluster potential 
Eq.(l) is fixed. One should note that it describes the interac- 
tion between only two clusters; if we consider other clusters' 
influences, a cluster will always be in the minimal potential 
state. When the cluster deviates from the equilibrium po- 
sition, it will be pulled back due to the stronger repulsion 
from one side, just as the case of a chain of springs. 

3.3 Comparison with the MIT bag model 

In the MIT bag model, quark matter is composed of mass- 
less up and down quarks, massive strange quarks, and few 
electrons. Quarks are combined together by an extra pres- 
sure, denoted by the bag constant B. For the comparison , 
we apply the formulae given by Alcock (|Alcock et al.lll986h 
to calculate the equation of state, with strange quark mass 
m s = lOOMeV, the strong coupling cons tant a B = 0.3 , and 
the bag constant B = 60MeV/fm -3 fe.g JZdunikkOOOh . The 
comparison of equation of state in our model and in the MIT 
bag model is shown in Fig.l. 

In our model, quarks are grouped in clusters and these 
clusters are non-relativistic particles. If the inter-cluster po- 
tential can be described as the Lennard-Jones form, the 
equation of state can be very stiff, because at a small inter- 
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Figure 1. The comparison of the equations of state, in the case 
N q = 3, including Uo = 50 MeV (blue solid lines) and Uo = 
100 MeV (blue dashed lines), and the corresponding case N q = 
18 with U = 50 MeV (red dash-dotted lines) and Uo = 100 
MeV (red dotted lines), and that derived in the MIT bag model 
with the mass of strange quark m B = 100 MeV and the strong 
coupling constant a B = 0.3 and bag constant B = 60MeV/fm -3 
(thin lines), for a given surface density p B = 2po- Here and in the 
following figures, po is the nuclear saturation density. 

cluster distance (i.e., the number density is large enough), 
there is a very strong repulsion. Whereas in MIT bag model 
quarks are relativistic particles (at least for up and down 
quarks). For a relativistic system, the pressure is propor- 
tional to the energy density, so it cannot have stiff equation 
of state. 

3.4 The speed of sound 

The adiabatic sound speed is defined as 

c s = ^dP/dp. (19) 

If we use the equation of state in our model, the speed of 
sound will exceed the speed of light not far away from the 
surface of a quark star. It seems to contradict to the rela- 
tivity that signals cannot propagate faster than light. 

The possibility of speed of sound exceeding the speed 
of ligh t in ultradense matter have been discussed previ- 
ously (jBludman fc Rudermanlll968l ) because of using clas- 
sical potential (i.e., a kind of action at a distance). The 
physical r easons of apparent superlu minal have also been 
analyzed (jCaporaso fc Brecherl Il979l ) . The authors argued 
that the adiabatic sound speed can exceed the speed of light, 
yet signals propagate at speed less than c. 

One reason is that the P(p) relation arises from a static 
calculation, ignoring the dynamics of the medium. The no- 
tion that c s is a signal propagation speed is a carry-over 
from Newtonian hydrodynamic, in which one assumes infi- 
nite interaction speed but finite temperature, so the static 
and dynamic calculations give the same result. On the other 
hand, if one assumes finite interaction speed and zero tem- 
perature, the adiabatic sound speed is not a dynamically 
meaningful speed, but only a measure of the local stiffness. 
Another reason is that a lattice does not have an infinte 
range of allowed frequencies of vibration, but a signal should 
contain components at all frequencies. Therefore, the adia- 



Figure 2. The mass-radius and mass-central density (rest-mass 
energy density) curves, in the case JVq = 3, including Uo = 50 
MeV (blue solid lines) and U = 100 MeV (blue dashed lines), 
and the corresponding case Nq = 18 with Uo = 50 MeV (red 
dash-dotted lines) and Uo = 100 MeV (red dotted lines), for a 
given surface density p B = 2po- 

batic sound speed is not capable of giving the velocity of 
propagation of disturbs. 

In our model, although we have not make it explicit 
that how the particles interact with each other, we may as- 
sume that the interaction is mediated by some particles with 
non-zero masses, and the interaction does not propagate in- 
stantaneously. We have also use the low frequency approx- 
imation to calculate the lattice energy. Therefore, we could 
conclude that in our model the signal can not propagate 
faster than light. 

Whether the equation of state of cold quark matter can 
be so stiff that the adiabatic speed of sound is larger than 
c could still be an open question. However, in our present 
paper we do not put the limitation on the adiabatic sound 
speed, and only treat it as a measurement of the stiffness of 
the equation of state. 



4 MASSES AND RADII 

From the equations of state, we can get the mass-radius 
curves and mass-central density curves (the central density 
only includes the rest mass energy density), as are shown in 
Fig.2. 

Because of stiffer equations of state, which we have dis- 
cussed in §3, the maximum masses of quark stars in our 
model could be higher. In Fig.2, we can see that (i) a deeper 
potential well Uo means a higher maximum mass; (ii) if there 
are more quarks in a quark-cluster, the maximum mass of a 
quark star will be lower. 

A stiffer equation of state leading to a higher max- 
imum mass could have very important astrophysical im- 
plications. Although we could still obtain high maximum 
mass es under MIT bag m odel by choosing suitable parame- 
ters <|Zdunik et alJ feoOO), a more realistic equati on of state 
in th e density-dependent quark mass model (e.g.. iDev et al.l 
|l998f) is very difficult to reach a high enough maximum stel- 
lar mass, which was considered as po ssible negative evidence 
for quark stars l|Cottam et al.l [20021 ). Some recent observa- 
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tions have indica ted some massive (~ 2Mq) pulsars (e.g., 
iFreire et alj|2008h ; however, because of the uncertain incli- 
nation of the binary systems, we are still not sure about the 
real mass. Though we have not definitely detected any pul- 
sar whose mass is higher than 2Mq up to now, the Lennard- 
Jones quark star model could be supported if massive pul- 
sars (> 2Mq) are discovered in the future. Moreover, a 
high maximum mass for quark stars might be helpful for 
us to understand the mas s-distribution of stellar-mass black 
holes (|Bailv n et al.|[l998l ). since a compact star with a high 
mass (e.g., ~ 5Mq) could still be stable in our model pre- 
sented. 



5 CONCLUSIONS AND DISCUSSIONS 

In cold quark matter at realistic baryon densities of com- 
pact stars (with an average value of ~ 2 — 3po), the interac- 
tion between quarks is so strong that they would condensate 
in position space to form quark-clusters. Like the classical 
solid, if the inter-cluster potential is deep enough to trap the 
clusters in the potential wells, the quark matter would crys- 
tallize and form solid quark stars. This picture of quark stars 
is different from the one in which quarks form cooper pairs 
and quark stars are consequently color super-conductive. 

In this paper, we argue that quarks in quark stars are 
grouped in clusters and the quark-clusters form simple-cubic 
structure, and apply Lennard-Jones potential to describe the 
interaction potential between quark-clusters. The parame- 
ters such as the depth of potential Uq (50 MeV and 100 
MeV) and the range of interaction ro (about 1 to 3 fm) are 
given by the physical context of quark stars. Under such 
equations of state, the masses and radii of quark stars are 
derived, and we find that the mass of a quark star can be 
higher than 2M Q . 

It is surely interesting to experimentally or observa- 
tionally distinguish between our solid quark star model 
and other models for quark stars, e.g., the color super- 
conductivity state. Starquakes could naturally occur in solid 
quark stars and the observations of pulsar glitches and SGR 
giant flares cou ld qualitatively be reproduced wh en the solid 
matter breaks l|Zhou et all |2004| ; IXu et alj |2006). moreover 
the post-glitch recovers in the solid quark star model and the 
color super-conductivity model would be different. Addition- 
ally, because the solid quark star model depends on quark 
clustering, the interaction behaviors between quarks could 
be tested in sQGP (strongly coupled quark-gluon plasma, 
see Shuryak 2009) by the LHC and/or FAIR experiments. 

The research of compact stars involves two kinds of 
challenges: particle physics and many-body physics. Nev- 
ertheless, if we know about the properties of compact stars 
from observations, we can get information of the elemen- 
tary physics. Take the model we discussed in this paper as 
an example. If we get the masses and radii of some pulsars 
from accurate enough observations, we can put limits on 
the parameters such as potential well depth Uo, interaction 
range ro, and the number of quarks that condensate in posi- 
tion space to form a cluster, which could help us to explore 
the strong interaction between quarks. Although the state 
of cold quark matter at a few nuclear densities is still an un- 
solved problem in the low-energy QCD, it would be hopeful 



for us to use pulsars as idea laboratories to study the nature 
of strong interaction. 

In general, stars are equilibrium bodies with pressure 
against gravity. The thermal and radiation pressure dom- 
inates in main sequent stars, while degenerate pressure of 
Fermions, originated from Pauli's principle, dominates in 
Fermion stars (e.g., while dwarfs). For solid quark stars in 
the models presented in this paper, the pressure is related 
to the increase of both potential and lattice vibration ener- 
gies as the stellar quark matter contracts. The degenerate 
pressure might be negligible there. 
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